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We present a new measurement of the mass of the top quark using lepton + jets tt events collected 
by the D0 experiment in Run I of the Fermilab Tevatron Collider. The mass is extracted through a 
comparison of each event with a leading-order matrix element that depends on the top quark mass. 
The result is Mt — 180.1 ± 3.6 (stat) ± 3.9 (sys) GeV/c^. Combining this improved measurement 
with our previous value from dilepton channels yields the new D0 result Mt — 179.0 ± 3.5 (stat) ± 
3.8 (sys) GeV/c^. 



The observation of the top (t) quark [J, ^ was one of the major confirmations of the validity of the standard 
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model (SM) of particle interactions. Through radiative 
corrections of the SM, the mass of the top quark (M^), 
along with that of the W boson (Mw) 13 , constrains the 
mass of the hypothesized Higgs boson Q . Mw is known 
to a precision of < 0.1%, while the uncertainty on Mt is 
at the 3% level 3] . Improvements in both measurements 
are required to further limit the mass range of the Higgs 
boson, and to check the self-consistency of the SM. It is 
therefore important to develop techniques for extracting 
a more precise value of M^. 

We report a new measurement of the mass of the top 
quark using tt events containing an isolated lepton and 
four jets, collected by the D0 experiment in Run I of 
the Fermilab Tevatron Collider. The data correspond to 
an integrated luminosity of 125 pb~^, and this analysis 
is based on the same data sample used to extract Mt in 
a previous publication |^. 

As before, we assume that the top quark decays 100% 
of the time to a W'^ boson and a b quark, which for a 
ti pair implies W'^W~'bb in the final state. This analysis 
is based on decay channels containing a lepton (electron 
or muon from one W — > Ivi decay ) and jets (from the 
evolution of the b quarks and the quarks from the other 
W — > qq' decay) in the final state p- After offline selec- 
tions on lepton transverse energy {e!^^ > 20 GeV) and 
pseudorapidities (It^^I < 1.7 for muons and \T]e\ < 2.0 for 
electrons), on jet transverse energies {Et > 15 GeV) and 
pseudorapidities (|?7| < 2.0), imbalance in transverse en- 
ergy (^T > 20 GeV), and on W boson decay products 
{El^^+f^T > 60 GeV) and pseudorapidity H-qw] < 2.0), 
and after applying several less-important criteria Q , the 
original event sample consists of 91 events with one iso- 
lated lepton and four or more jets. (Unlike the previous 
analysis, we do not distinguish between events that have 
or lack a muon associated with one of the jets, signi- 
fying the possible presence of a 6-quark jet in the final 
state.) The new analysis involves a comparison of these 
91 events with a leading-order matrix element for tt pro- 
duction and decay. To minimize the effect of higher-order 
corrections, we restrict the study to events containing ex- 
actly four jets, which reduces the sample to 71 events. 

In the previous analysis, the four jets with highest 
Et were assumed to represent the four quarks in the 
event. These, along with the lepton and the unob- 
served neutrino were fitted to the kinematic hypothesis 
pp ^ ti ^ WWbb, subject to the constraints of overall 
momentum-energy conservation, the known mass of the 
W boson, and the fact that the unknown mass of the 
top quark was assumed to be identical for the top and 
antitop quarks in the event. With twelve ways to per- 
mute the jets, there were twelve possible fits (six when 
one of the jets was tagged as a 5 jet), and the solution 
with lowest was chosen as the best hypothesis, thereby 
defining the fitted mass rngt for the event. The same 
procedure was used to generate templates in variables 
of interest as a function of input top quark mass. This 
was based on the HERWIG Monte Carlo (MC) program 
0, which was used to generate events that were passed 



through full detector simulation and event reconstruction 
l8j|. Background events, consisting mainly of multijets 
(20%) and W^-|-jets (80%) production, were processed in 
a similar manner. The background from multijet pro- 
duction was based on studies of multijet events in data 
0, and the background from M^-l-jets events was based 
on events generated with VECBOS A four- variable 
discriminant (D) defined the probability that an event 
represented signal as opposed to background. A proba- 
bility density was defined as a function of the discrim- 
inant D and mgt, and a comparison of data and MC 
via a likelihood was used to determine the most likely 
mass of the top quark. The resulting measurement is 
Mt = 173.3 ± 5.6 (stat) ± 5.5 (sys) GeV/c^. 

The new method is similar to that suggested for tt 
dilepton decay channels [loj | , and used in a previous mass 
analyses of dilepton events [ll|. A similar approach has 
also been suggested for the measurement of the mass of 
the W boson at LEP 12]. Given TV events, the top quark 
mass is estimated by maximizing the likelihood: 

N 
i=l 

where Xi is a set of variables needed to specify the zth 
measured event, Pm is the probability density for observ- 
ing that event, and a represents the parameters to be de- 
termined (in this case a is the mass of the top quark) . De- 
tector and reconstruction effects are taken into account in 
two ways. Geometric acceptance, trigger efficiencies, and 
event selection enter through a multiplicative function 
A{x) that is independent of a, and relates the observed 
probability density Pm{x, a) to the production probabil- 
ity P(a;, a): Pm{x,a) = A{x)P{x^a). Energy resolution 
and merging and splitting of jets are taken into account 
in a "transfer" function, W{y,x), discussed below. The 
production probability density can be written as a con- 
volution of the calculable cross section and W^(?/,x): 

P{x,a) = [ da{y,a)dqidq2f{qi)f{q2)W{y,x) 

cr(a) J 

(2) 

where W{y,x), our general transfer function, is the nor- 
malized probability density that the measured set of vari- 
ables X arise from a set of partonic variables j/, da{y, a) 
is the partonic differential cross section, and f{qi) are 
parton distribution functions for the incoming partons 
with longitudinal momenta qi. Dividing by cr(a), the to- 
tal cross section for the process, ensures that P{x, a) is 
properly normalized. The integral in Eq.ljSJ sums over 
all possible parton states leading to what is observed in 
the detector. 

For the tt production probability, the measured angles 
of the jets and of the charged lepton are assumed to be 
the angles of the partons in the final state. Given the 
detector resolutions the electron energy is assumed to be 
exact, and the muon energy is described by its known 
resolution [l3. Evaluation of Eq.® for the e-|-jets chan- 
nel involves two incident parton energies (we take these 
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partons to be quarks, and ignore the «10% contribution 
from gluon fusion), and six objects in the final state. 
The integrations over the essentially fifteen sharp vari- 
ables (three components of electron momentum, eight jet 
angles, and four equations of energy-momentum conser- 
vation), leave five integrals that must be performed to 
obtain the probability that any event represents tt pro- 
duction for some specified value of top quark mass Alt: 



E 

pcrm,i/ 



dpidmldMldmldM"^ 
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For [A^tfP, we use the leading-order matrix element f^li 
f{qi) and /(g2) are CTEQ4M parton distribution func- 
tions for the incident quarks $6 is the phase-space 
factor for the six-object final state, and the sum is over 
all twelve permutations of the jets (the permutation of 
the jets from W boson decay was performed by sym- 
metrizing the matrix element), and the up-to-eight pos- 
sible neutrino solutions. Conservation of transverse mo- 
mentum is used to calculate the transverse momentum 
of the neutrino. Wjcts(£^part, -Ejet) is the part of W{y,x) 
that refers to the mapping between parton-level ener- 
gies impart and energies measured in the detector -Ejct- 
Four of the variables chosen for integration (mi. Mi, 
7712 and M2), namely the masses of the W bosons and 
of the top quarks in the event, are economical in com- 
puting time, because the value of |A^ttP is essentially 
negligible except at the peaks of the four Breit-Wigner 
terms in the matrix element, pi is the energy of one 
of the quarks in the hadronic decay of one of the W 
bosons. The narrow-width approximation is used to in- 
tegrate over the top quark masses, and Gaussian adaptive 
quadrature is used to perform the three remaining 
integrals. VFjots(^'part, ^jot) is the product of four fimc- 
tions i^(i?pai.t J -Ejot)' '^^^ for each jet, with a functional 
form of the sum of two Gaussians, with parameters hav- 
ing linear dependence on -Ep^it • The parameters used for 
6 quarks are different from those for the lighter quarks, 
and there are therefore twenty jet energy parameters in 
all. About 15,000 simulated tt events (generated with 
masses between 140 and 200 GeV/c^ in HERWIG, and 
processed through detector simulation) are used to deter- 
mine the above twenty parameters. For a final state with 
a muon, VFjets is expanded to include the muon momen- 
tum resolution, and an integration over muon momentum 
is included in Eq.©. 

The W+A jets matrix element from VECBOS is used 
in Eq.Q to calculate the background probability Pbkg- 
The integration is performed over the energy of the four 
partons leading to jets and the M^-boson mass. The prob- 
ability is summed over the twenty-four jet permutations 
and two neutrino solutions. The integration over parton 
energies is performed using MC techniques, increasing 
the number of random points until the integral converges. 
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Figure 1: (a) Distribution in probability of events being back- 
ground, and (b) discriminant Pu / (Ptt + Pbkg) , calculated for 
the 71 tt candidates (data points). The data are compared 
with results expected for the sum (open histogram) from 
MC-simuIated sources of tt (left-hatched) and lF-|-jets (right- 
hatched) events. Only events with Pbkg < 10~^^ (indicated 
by the vertical line) are considered in the final analysis. 



(MC studies show that the 20% background from multijet 
events is represented satisfactorily by that for VF-|-jets.) 

After adding the probabilitites for the non-interfering 
tt and W+A jets channels, the final likelihood as a func- 
tion of Mt is written as: 

JV 

- lTLL{a) = - ^ ln[ciPtiixi,a) + C2Phkg{xi)] 

i=l 

+ Nci / A{x)Pu{x, a)dx + Nc2 / A{x)P^kg{x)dx 



The above integrals are calculated using MC methods, 
for which the acceptance A{x) is 1.0 or 0.0, depending 
on whether the event is accepted or rejected by the anal- 
ysis criteria. The best values of a, representing the most 
likely Mt , and the parameters Ci are defined by minimiz- 
ing — In L(q). 

Studies of samples of HERWIG MC events used in the 
previous analysis indicate that the new method should 
yield almost a factor of two reduction in the statistical 
uncertainty on the extracted Mt- These studies also re- 
veal that there is a systematic shift in Alt that depends on 
the amount of background in the data sample. For high 
statistics, the shift is about 2 GeV/c^ when the back- 
ground approaches 80% of total. To minimize this bias, 
a selection is introduced based on the probability that an 
event represents background from M^-|-jets. Figure ^a) 
shows a comparison between the probability for a back- 
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Figure 2: (a) Negative of the log of the likelihood as a function 
of the top quark mass, (b) The likelihood normalized to its 
maximum value in plot (a). The curve is a Gaussian fit to the 
likelihood plot. The hatched area corresponds to the 68.27% 
probability interval. 



Figure 3: Mass of the top quark as a function of the cutoff in 
background probability. The number of remaining events is 
shown above each point. The point with the larger dot is the 
value used in this analysis. 



ground interpretation of a large sample of mixed MC 
events (upper-most histogram) and the 71 tt candidates 
(data points). The total number of MC events is nor- 
malized to the 71 4-jet ti candidates. The left-hatched 
(right-hatched) histogram shows the contribution from tt 
{W + A jets) MC events to the total. Only the 22 events 
to the left of the vertical line are chosen for the final 
analysis (Pbkg < 10~"). The ratio oi tt to W + 4 jets 
events in the MC is normalized to the 12/10 ratio found 
for the data to the left of the vertical line, as described 
below. (The selected value of Pbkg < 10~^^ is based on 
MC studies carried out before applying the method to 
data, and, for a top quark mass of 175 GeV/c^, it retains 
71% of the signal and 30% of the background.) 

A discriminant D = P^f/^P^^ + Pbkg) was defined to 
quantify the likelihood for an event to correspond to sig- 
nal ^. Since the signal probability depends on Mt, D 
was calculated with the signal probability taken at its 
most likely value. Figure^b) shows a comparison of the 
discriminant calculated for data (points with error bars) 
and for MC events (open histogram), with the MC nor- 
malized as in Fig. ^^a). Since the discriminant depends 
directly on Mt, it was not used to reject background and 
is shown simply to illustrate the level of discrimination 
of signal from background. 

FigureEta) shows the value of — InL(a) as a function 
of Mt for the 22 events that passed all selection criteria. 
— InL(a) was minimized with respect to the parameters 
Ci at each mass point. Figure |2b) shows the likelihood 
normalized to its maximum value. The Gaussian fit in 
the figure yields Mt = 179.6 GeV/c^, with an uncer- 
tainty 5Mt = 3.6 GeV/c^. MC studies show that [l7j : 
(i) SMt is compatible with the uncertainties obtained 
in MC ensemble tests, and (ii) there is a shift of -0.5 
GeV/c^ in the extracted mass. After applying the 0.5 
GeV/c^ correction, our new value of the top quark mass 
is Mt = 180.1±3.6 (stat) GeV/c^ As Fig. [Ja) indicates, 
the cutoff chosen in Pbkg does not reduce significantly the 
number of tt events, and therefore Mt should be stable 
relative to variations in this cutoff. Figure |31 shows that 



a change in the cutoff in Pbkg by more than an order 
of magnitude changes the number of events used in the 
analysis by more than a factor of two, but, as expected, 
does not have a significant impact on Mt. 

The total number of tt events in the 91-event sample 
is deduced to be (11 ± 3)/(0.71 x 0.70 x 0.87) = 25 ± 7, 
where 11 is the number of extracted (using ci and C2) 
events, and the corrections are for: (?) acceptance (0.71), 
(m) events with more than 4 jets (0.70), and {in) tt 4- 
jets events that appear as background because a leading- 
order matrix element does not represent them correctly 
(0.87). The number of tt events is consistent with that 
found in the previous analysis j^. The uncertainties in 
the quoted efficiencies are negligible compared with the 
statistical uncertainty in the total number of tt events. 

In the previous analysis 0, 7-l-jet events were used 
to check the energy scale in the experiment relative to 
MC simulation. This calibration had an uncertainty of 
SE= 0.025 E + 0.5 GeV. Consequently, we rescaled the 
energies of all jets in our sample by ±SE, redid the anal- 
ysis, and the average of the two rescaled results for Mt 
{SMt = (3.0 + 3.5)/2 « 3.3 GeV/c^) is taken as the sys- 
tematic error due to the uncertainty in the jet energy 
scale (JES). Additional contributions to the systematic 
uncertainty are listed in Table HI and details on the eval- 
uation of these systematic errors and further description 
of the analysis technique can be found in Ref. 17]. 

Our method provides substantial improvements in 
both the statistical and systematic uncertainties. This 
is due to two main differences relative to the previous 
analysis: («) each event now has an individual probabil- 
ity as a function of the mass parameter, and therefore 
well-measured events having a narrower likelihhod con- 
tribute more to the extraction of the top quark mass than 
those that are poorly measured, and {ii) all possible jet 
and neutrino combinations are included, which guaran- 
tees that all signal events contribute to the measurement. 

In conclusion, we have presented a new measurement of 
the mass of the top quark using a method that compares 
each individual event with the expected differential cross 
section for tt production and decay. We obtain a sig- 
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Table I: Systematic uncertainties on the measurement of Mt- 



iviociei ror ll 


1.1 LtC V / C 


Model for backgound(14^-t-jets) 


l.U LtC V / C 


Noise and multiple interactions 


i.o LreV 1 C 


Jet energy scale 


3.3 GeV/c^ 


Parton distribution function 


0.1 GeV/c^ 


Acceptance correction 


0.3 GeV/c^ 


Bias correction 


0.5 GeV/c^ 


Total 


3.9 GeV/c^ 



nificant improvement in statistical uncertainty over the 
previous measurement that is equivalent to a factor 
of 2.4 more data. 

From the differences in the two analyses, and from sta- 
tistical fluctuations arising from using a subsample of the 
original data, we expect the difference between the orig- 
inal and the new mass measurement to be on the order 
of 4 GeV/c^. Thus, the two results differ by less than 
two standard deviations. The current analysis is also 
less sensitive to the calibration of the JES, and leads to 
an improved systematic uncertainty. The new result is: 

Mt = 180.1 ± 3.6 (stat) ± 3.9 (sys) GeV/c^. 

Combining the two uncertainties in quadrature, we ob- 
tain Ml = 180.1 ±5.3 GeV/c^, which has an uncertainty 
comparable to all the previous measurements of D0 and 



CDF Hi combined. 

Using the procedure described in Ref.fTsIl, the new 
measurement can be combined with that obtained us- 
ing the dilepton sample collected at D0 during Run I 
yielding the new D0 value for the mass of the top 
quark: 

Mt = 179.0 ± 3.5 (stat) ± 3.8 (sys) GeV/c^ 

This is the most accurate measurement of the top quark 
mass in any single experiment. The impact of the new 
D0 top-quark mass measurement on the world average 
top-quark mass as well as on Higgs and supersymmetry 
constraints is a subject of a separate recent publication 
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